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Abstract. - Many microorganisms swim through gels, materials with nonzero zero-frequency 
elastic shear modulus, such as mucus. Biological gels are typically heterogeneous, containing 
both a structural scaffold (network) and a fluid solvent. We analyze the swimming of an infinite 
sheet undergoing transverse traveling wave deformations in the "two-fluid" model of a gel, which 
treats the network and solvent as two coupled elastic and viscous continuum phases. We show 
that geometric nonlinearities must be incorporated to obtain physically meaningful results. We 
identify a transition between regimes where the network deforms to follow solvent flows and where 
the network is stationary. Swimming speeds can be enhanced relative to Newtonian fluids when 
the network is stationary. Compressibility effects can also enhance swimming velocities. Finally, 
microscopic details of sheet-network interactions influence the boundary conditions between the 
sheet and network. The nature of these boundary conditions significantly impacts swimming 
speeds. 



While biological locomotion at low-Reynolds number is 
a well-established and active field (see [T] for a review), 
only recently has swimming in complex, non-Newtonian 
media started to be systematically explored. Indeed, many 
microorganisms routinely swim through complex or non- 
Newtonian media. Bacteria in the stomach such as He- 
licobacter pylori encounter gastric mucus [2] ; mammalian 
sperm swim through viscoelastic mucus in the female re- 
productive tract [Sill]; and spirochetes burrow into the 
tissues they infect [5]. Most of the prior theoretical work 
on swimming in complex media has focused on viscoelastic 
fluids [6Ul4]. Because prescribed swimming strokes lead 
to the same swimming speed in linearly viscoelastic fluids 
as in Newtonian fluids [ZlIIS], these studies have focused 
on the effects of nonlinear viscoelasticity [6ll8ll9llTTHl4] and 
the effect of altered medium response on swimming strokes 
via fluid-structure interactions [9lll0j. However, many bi- 
ological viscoelastic materials contain crosslinked polymer 
networks, and are therefore solid rather than fluid. In 
this paper we focus on swimming in such crosslinked ma- 
terials, which have a nonzero zero-frequency elastic shear 
modulus. We use the term "gel" to refer to these types of 
materials. 

Generally, viscoelastic materials have frequency- 
dependent response, while the distinction we make be- 
tween solid and fluid is a zero- frequency property. Some of 



the issues addressed by studies of viscoelastic fluids carry 
over straightforwardly to gels. For example, the fluid- 
structure interactions which determine the beating pat- 
terns of sperm flagella depend on the viscoelastic response 
at the beating frequency, not at zero frequency. However, 
a simple physical argument demonstrates a fundamental 
difference between swimming in a gel versus swimming in 
a fluid. In fluid dynamics, one uses no-slip boundary con- 
ditions, but for a swimmer in a solid, no-slip boundary 
conditions do not allow swimming. Under no-slip bound- 
ary conditions a swimmer would drag solid material along 
with its motion, leading to large deformations and restor- 
ing forces which oppose translation. Thus, swimming in a 
gel requires that the swimmer slips past the medium. This 
simple consideration highlights the importance of bound- 
ary conditions for understanding swimming in gels, a ma- 
jor theme of this paper. 

Another motivation for studying swimming in gels is 
the search for mechanisms that enhance swimming speeds 
in complex materials. Due to the observation that spiro- 
chetes can swim faster in fluids with an elastic response 
component [15j . there has been an interest in how vis- 
coelastic media can increase or decrease swimming ve- 
locities. For filamentous geometries, it has been sug- 
gested that polymeric fluids may have higher effective 
anisotropics in resistive-force theories, leading to increased 
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swimming speeds |16lll7j . Furthermore, the breaking of 
kinematic reversibihty (TB] in nonhnearly viscoelastic flu- 
ids makes certain reciprocal swimming strokes more ef- 
fective than in Newtonian fluids [TT1[T3]. Restricting at- 
tention to sheet-like geometries undergoing traveling- wave 
deformations (the subject of this paper), an infinite sheet 
in a nonhnearly viscoelastic fluids swims more slowly than 
in a Newtonian medium [8]. However, finite- length sheet- 
like swimmers in viscoelastic fluids I14| , and infinite swim- 
mers in a Brinkmann fluid can show increased swimming 
velocities [TT l fTQ ] . 

The latter two studies [TTlfT^ are closest in spirit to this 
letter. Both investigate how Taylor's swimming sheet [20] . 
an infinite sheet with a small amplitude traveling wave de- 
formation, behaves in the Brinkmann model, which treats 
the effect of a stationary porous phase on fluid flows in 
an effective medium approach. Instead, we solve for the 
flows induced by the moving sheet in the "two-fluid model" 
description of a gel [2TI - I24] . which represents the gel as 
two phases. One phase describes an elastic polymer net- 
work with displacement field u and Lame coefhcients ^ 
and A, and the other phase describes a viscous fluid sol- 
vent with velocity field v and viscosity rj. The two phases 
are coupled by a drag force proportional to their relative 
local velocity, with drag coefficient F, giving the governing 
equations 
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Here the network and solvent stress tensors are <t" and cr^, 
respectively. The moduli in the network stress tensor are 
effective, macroscopic moduli ~ for example, they incor- 
porate osmotic effects, and they depend implicitly on the 
volume fraction of the network. Eq. [5] and [H] express vol- 
ume conservation for each phase in terms of the network 
volume fraction ip assuming that the individual solvent 
and network phases are incompressible. Combining the 
two yields 



(7) 



In the rest of the paper, we work in the regime where 
network volume fraction (p can be ignored, but yet the 
macroscopic network stresses cannot be ignored, i.e. if 
we expand the solutions in ip, we obtain the zeroth-order 
solutions in tp. In this tp ^ limit, Eq. [7] reduces to 
V • V = 0. The drag coefficient F can be estimated as 
F « v/S.'^j where ^ is the mesh size of the network [21]. 




Fig. 1: Effective sheet-network boundary conditions depend on 
the microscopic swimmer-network interactions, a) Network is 
dense on the scale of swimmer deformations leading to direct 
network-sheet boundary conditions (Eqs. [8] and [Qjl. b) Net- 
work is sparse on the scale of swimmer deformations leading to 
solvent-mediated network-sheet boundary conditions (Eg. llOp . 



Finally, note that the network stress tensor incorporates 
nonlinear terms (involving products of u) explicitly in d^". 

Our use of the two-fluid model explicitly recognizes the 
role of heterogeneity in complex media, while still retain- 
ing a continuum description. Ref. |17j emphasizes that 
polymers in a viscoelastic fluid may be viewed as heteroge- 
neous inclusions in a viscous fluid. Treating the inclusions 
as a stationary background results in the Brinkmann fluid 
as an effective model. The two-fluid model allows us to 
describe the heterogeneous phase more realistically as a 
dynamic component, and understand when it is appropri- 
ate to think of heterogeneities as stationary or dynamic. 

Our investigation leads to three main conclusions. First, 
network dynamics determine when the presence of hetero- 
geneities leads to increased swimming velocities. Roughly 
speaking, network deformations are a result of drag forces 
arising from solvent flow, and are opposed by the stiffness 
of the network. For fast flows or compliant networks, the 
networks move with the flows and do not attenuate the 
flow, while for slow flows or stiff networks, the network 
acts like a stationary phase, leading to enhancements of 
the swimming velocity. We also find that compressibil- 
ity of the network can lead to additional enhancement of 
swimming velocities. 

Second, a consistent treatment of swimming velocities 
requires proper treatment of nonlinearities. We show that 
nonlinearities in the constitutive relationship (o-" in Eq. [3]) 
of the elastic phase of a gel do not affect swimming speeds. 
However, geometric nonlinearities arise from the terms of 
Eqs. [T] and [5] that involve the velocity of material points 
of the network, ^u. The displacement field is a function 
of current coordinates - i.e., the position of a material 
element originally at xq is x(xo,t) = xq + u(x, i) - so 
the velocity of material elements is given by the material 
derivative, defined implicitly via = dtu+ • Vu. In 
the remainder of the paper we refer to the term • Vu 
as the "convective" nonlinearity. Without it one cannot 
obtain physically meaningful swimming velocities. 

Finally, we highlight the importance of boundary condi- 
tions (BC) between the sheet and elastic phase. There is a 
long history of experimental support for the validity of no- 
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slip boundary conditions for fluids, and we use no-slip BC 
between the sheet and solvent, v(x, h{x, t), t) = dth{x, t)y, 
where h{x, t) is height of the sheet at (x, t). However, there 
is no a priori selection of boundary conditions between the 
sheet and network in our macroscopic theory. In this pa- 
per, wc consider two types of network-sheet BC that may 
plausibly result from microscopic interactions at a swim- 
mer surface. 

The first class of boundary conditions applies when the 
network is dense on the scale of the swimmer (Fig. [T^). 
Then we imagine that the sheet presses directly against 
the polymer network and the network must conform to the 
sheet. Using reference positions for the network produced 
by a flat sheet, 

Uy{x,h{x,t),t) ~ h{x,t). (8) 

The remaining boundary condition is applied to the tan- 
gential network stress via a kinetic friction law, with tan- 
gential stress proportional to the tangential velocity of the 
polymer relative to the sheet: 

t • (T'\x,h{x,t),t) ■ n = 

St • (^^u{x, h{x, t), t) - dth{x, t)y^ (9) 

where t and n are the tangent and normal to the sheet. 
Note that S — >■ corresponds to free tangential slip, while 
S — cxD corresponds to no-slip. 

In the second class of boundary condition, the network 
is dilute on the scale of the swimmer (Fig. [TJj) , so the sheet 
exerts no forces directly on the network. Instead all forces 
on the network are mediated by the solvent. Then the two 
components of the network stress acting on the sheet must 
vanish: 

a''{x,h{x,t),t) -h^O. (10) 

Both the "direct" BC (Eqs. [8l9| and "solvent-mediated" 
BC (Eq. [T0| allow the swimmer to slide past the network, 
in accordance with our physical discussion above. 

Solutions. In the frame where the sheet is stationary, 
the transverse displacements are described by h(x,t) = 
Re[&exp(i(A;a; — ^t)]- Throughout this paper we will be 
concerned with the zero-Reynold's number limit appro- 
priate for swimming microorganisms, and hence always 
ignore inertial effects. 

The sheet swims with velocity V,, but sufficiently far 
from the sheet in the lab frame the elastic phase is sta- 
tionary. Thus we write the displacement field in the frame 
of the sheet as u(x, t) ~ Vstx + 6u{x,t), with 6u small. 
Note that Vs turns out to be second order in the displace- 
ments. The solution strategy follows Taylor [20]. For small 
traveling wave amplitudes relative to wavelengths, we de- 
termine the velocity and displacement fields order by order 
in b. We write (for example) 6u = Ju*-^-* + Su'^^^ + with 
5u^^^ first order in b, etc. The fields satisfy the governing 
equations and boundary conditions, also expanded order 
by order in b. 



The first-order governing equations are 
^V2,5u(i) + (A + A*)V(V- W^)) = r(5t(5u(i) - v^i)) 

The first-order solutions take the form v'^^ ~ 
Re{v(y) exp(i(fca; — ujt))} (for example). The first order 
solutions can be found by introducing stream and poten- 
tial functions, v'^^ = curl(5'z) and Su^^^ = curl($z) -I- Vx, 
and then following the methods detailed in . Once the 
first order solutions are obtained, the swimming velocity 
is found by solving the second order equations. The swim- 
ming velocity is determined by x-component of the time- 
and x-averaged velocity and displacement fields. Note 
that any overall time- or x-derivative vanishes upon av- 
eraging, so the time-averaged z-componcnts of the 2nd 
order pieces of Eq. [1] and [2] are 

= r[F.-(42)) + ((a,<5uW.v)<5uW)]. (11) 

Here the brackets denote time-and x-averaging and ay^^ 
includes the 2nd order contributions to the network stress 
involving products of two first order fields. The solu- 
tions to these second order equations are (wi ) = Vs + 

ci cxp{- ^/TJ^y) + f{y) and (Sui^^) = -(r?/M)(wi^^) +^(2/), 
where f{y) is the inhomogcneous solution satisfying [rjdy — 

T)f = -F(at(5u(i) • ViSwi^') and g{y) is the inhomo geneous 
solution satisfying fidyg = — (tr^y^''). Boundary conditions 
at infinity rule out other possible solutions. 

These solutions must satisfy the second order boundary 
conditions, obtained by expanding Eqs. [5] and [TU] using 
h y — dxhic and t x + dxhy. For direct sheet-network 
interactions, the relevant BC are 

see Eqs. [HI and [HI 

with all quantities evaluated at y = 0. Eqs. [12] and [13] 
also apply for solvent-mediated BC, provided that we set 
S = in Eq. [13] In Eq. [13] note that all pieces involv- 
ing i^xl^^) cancel due to the form of {5u^^). Therefore 
nonlincaritics in the constitutive relation for the polymer 
stress (i.e., o-") do not affect the swimming velocity at sec- 
ond order. Similar cancellations help to show that these 
solutions satisfy overall force balance on the swimmer. 

Eliminating ci from Eqs. 1121 and 1131 vields Vs in terms 
of first order quantities. For direct network-sheet BC, 

see Eq. [14] 

The fourth term on the right hand side, {dtSu^^'' ■VSul^''), 
arises from the convective nonlinearity. For the solvent- 
mediated BC, the swimming velocity is given by Eq. [14] 
with S = 0. However, the swimming velocity with solvent- 
mediated BC is not just the swimming velocity with direct 
BC for S = 0, since the first-order solutions differ. 
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(dthd^h) + {hdydtSuJ^'') + {dtSuy^^^d^h) + (at^u^i) • ysu^^^ 
^ {{hdyv^^^) + /(o)) + |a,/(o) - h {{hdyu'^Ji^) + {dM^l'i^ - ^Si'^)> 



(14) 



For an incompressible network (A — > oo) the swimming 
velocity has a simple form when the convectivc nonlinear- 
ity is included: 



3/2 
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rr/[4 + 6^r/(7^fc2) + r/(7?fc2)] 

[E + ^)[2 + .jTjW)Y 



(15) 



For an incompressible network the swimming velocity is 
the same as in a Newtonian fluid, Vn ~ 6^a;fc/2, when the 
convective nonlinearity is not included. 

Results and discussion. Our solutions reveal that 
the heterogeneous nature of gels leads to behavior not 
captured by modeling viscoelastic media as single-phase 
fluids. This is apparent even in the first order solutions. 
In a viscoelastic fluid, the first order solutions are the same 
as in a Newtonian fluid, while in our solutions, once the 
network becomes compressible, or once one uses solvent- 
mediated sheet-network BC, there is relative motion be- 
tween the solvent and network and both phases move dif- 
ferently from a Newtonian fluid. 

Turning to the swimming velocities, we see the impor- 
tance of correctly accounting for nonlinearities. As in pre- 
vious studies of the role of nonlinearities in single-phase 
viscoelastic fluids [8l[9l[Tl], for small amplitude swimming 
strokes the swimming speed is a second order effect, and so 
a consistent description must take into account any non- 
linearities in constitutive laws. In the case of a gel de- 
scribed using a two-fluid model, we have shown that for 
both direct and solvent-mediated sheet-network BC, non- 
linearities in the constitutive law for the elastic phase do 
not affect the swimming velocity. However, the geomet- 
ric convective nonlinearity must be taken into account to 
obtain physically meaningful swimming velocities. Fig. [2] 
shows the swimming velocity as a function of inverse fric- 
tion 1/S. Previously we argued that no-slip sheet-network 
BC prevent a sheet from swimming, which implies that the 
swimming velocity should vanish as l/S — > 0, (the no-slip 
limit). The figure shows that this expected behavior is 
only obtained when the convective nonlinearities are in- 
cluded. This holds for incompressible networks (A — > oo) 
as well as stiff (small ijlo/ ^l) and compliant (large rju)/ ^) 
compressible networks. In the rest of this paper we always 
include the convective nonlinearity. 
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Fig. 2: Swimming speed normalized by Newtonian swimming 
speed, Vs/Vn, as a function of rfk/'E. Only when convective 
nonlinearities are included does the swimming velocity van- 
ish for no-slip BC {rjk/'E. — >■ 0). All curves are for direct 
network-sheet BC. Thin-black (dots) and thick-black (solid) 
incompressible network with V /[-qk^) = 1. Red (long dashes) 
r?cj//i = 0.01, A = ^, T/l'qk^) = 10. Magenta (dash-dots) 
r\ujl\i = 10, A = /i, r/(77fc^) — 10. Blue (short dashes) 
rfu)/^ = 0.01, \ — fi, r/(?yfc^) = 10. Green (medium dashes) 
■quj/^i = 10, \ = ^i, r/(77fc^) = 10. Thick-black, Red, and 
Green curves include the convective nonhnearity; thin-black, 
magenta, and blue curves do not. 



Compressibility effects are a possible mechanism for 
swimming speed enhancement. Fig. |3] show the swim- 
ming speed as a function of Lame coefficient A, a mea- 
sure of ineompressibility (the bulk modulus of compression 
is 2fj, + A). For direct sheet-network BC, the strongest 
enhancement of swimming speed occurs for a material 
with large shear modulus (/i — > oo) but which is max- 
imally compressible (A = 0), and for a network which 
slips freely past the swimming sheet (S = 0). Compress- 
ibility does not strongly affect the swimming speed for 
solvent-mediated sheet-network BC. Swimming speeds can 
be faster than in Newtonian fluids for stiff networks/low 
frequencies (small rjuj/^), but not for compliant net- 
works/high frequencies (large rjuj/fj,), reflecting a broader 
distinction which we address next. 

Fig. |4] plots the swimming velocity as a function of 
?/a;//i. There is a transition between two regimes, one for 
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Fig. 3: Swimming speed normalized by Newtonian swimming 
speed, Vs/Vn, as a function of incompressibility, X/ fJ.. Com- 
pressibility effects can enhance swimming speeds for direct 
network-sheet BC. Swimming speed is enhanced relative to 
the Newtonian case for stiff networks or low frequencies (small 
riuj/y). Red and green curves are for direct network-sheet BC 
and H = 0; blue and magenta curves are for solvent-mediated 
network-sheet BC. Red (long dashes): r]U)/ ^ = 0.01, V/ij^k^) = 



Fig. 4: Swimming speed normalized by Newtonian swimming 
speed, Vs/Vn, as a function of rjuj/ ^. There is a transition 
at r]uj/fi ~ between fast swimming (stiff network/low 

frequencies, small r]u/ii) and slow swimming (compliant net- 
work/high frequences, large rju/fi). Black, red and magenta 
curves are for direct network-sheet BC and H = 0; blue and 
green curves are for solvent-mediated network-sheet BC. Black 
(solid): \ = ^i, T/{r]k^) = 10. Red (long dashes): A = 0, 



10. Green (medium dashes): rJUJ/^I = 10, r/(r?fc^) = 10. Blue r/{rjk^) = 10. Magenta (dash-dots): A = 0, r/{rjk^) = 1. 



(short dashes): rju/fi = 0.01, r/{r]k^) = 10. Magenta (dash- 
dots): Tjuj/fj. = 10, r/{rjk^) = 10. 
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flexible networks or high frequency beating (small rjuj / fi) , 
leading to a smaller swimming velocity; and one for rigid 
networks or low frequency beating (large rjuj/fi), leading 
to larger swimming velocities. The swimming speed en- 
hancement is larger for larger values of T, and is larger 
for solvent-mediated network-sheet BC than for direct 
network-sheet BC. Physically, the transition between the 
two regimes corresponds to whether the network does not 
deform (stiff network/low frequencies) or whether the net- 
work deforms to move together with the solvent (compli- 
ant network/high frequencies). Since the network deforms 
due to drag forces from the solvent, one can predict that 
the transition between the two regimes occurs when the 
force/volume due to drag forces, wuT, becomes compa- 
rable to the force/volume due to elasticity, ^ jiuk^, i.e. 
when rjui/ fi ~ ijk"^ /T, which can be confirmed by compar- 
ing the curves in Fig. |4l 

We can compare our results to those for the single-phase 
Brinkmann fluid [TTllin], which obeys rjS/^v — Vp = F^v 
and has swimming speed Vb = Vat-^/I + F-^/ {rjk'^). In 
the Brinkmann fluid the elastic phase is not a dynamic 
degree of freedom, but instead enters as an inert back- 
ground retarding solvent motion. One expects that the 
two-fluid model behaves like a Brinkmann fluid when the 
network does not deform (i.e., rjuj/fi ^ rjk^ /T) and if 
F^ ~ F. This expectation is borne out in Fig. [5l but 
only for solvent-mediated BC. When the network deforms 
(compliant networks/high frequencies) the network-sheet 
BC do not greatly affect the swimming speed (magenta 
and green curves). 

Our results highlight the importance of understanding 



the proper boundary conditions for swimmers in gels. The 
boundary conditions need not be the same for all swim- 
mers; they may depend on factors including swimmer mor- 
phology, surface biochemistry, or gel structure. However, 
one way to estimate whether a swimmer has direct inter- 
action with the polymer network in a gel is to compare the 
swimmer size with the mesh size of the network. For ex- 
ample, sperm swimming in cervical mucus have head sizes 
of 3-5 fim. The mesh size of the network of mucin fibers in 
cervical mucus varies strongly with time during the men- 
strual cycle, with mesh sizes of up to 25 ^m around ovula- 
tion, and as small as 200 nm at other times p5ll26| . Thus, 
during ovulation, when the sperm may be able to swim 
through the mucus while rarely encountering mucin fibers, 
solvent-mediated network-sheet BC likely apply. When 
the mesh sizes are small, direct network-sheet BC are more 
appropriate. Therefore, both types of boundary condi- 
tions described in this work may be relevant in biological 
situations. One important lesson is that even for an effec- 
tive medium model of a heterogeneous material (such as 
our two-fluid model), microscopic heterogeneity will have 
strong consequences for the boundary conditions and cor- 
responding swimming properties, since different bound- 
ary conditions lead to nearly 100% changes in swimming 
speeds in the stiff network limit. 

Note that for the particular case of sperm in cervical 
mucus, the sperm are likely to be in the compliant net- 
work/high frequency phase: for well-hydrated mucus (near 
ovulation) when the mesh size is « 10 //m, the shear mod- 
ulus is « 1 Pa, and viscosity is ~ 0.4 Pa s; while far- 
ther from ovulation with mesh size « 0.5 /xm, the shear 
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Fig. 5: Swimming speed normalized by Newtonian swimming 
speed, Vs/Vjv, as a function of drag coupling F/rjk. The swim- 
ming speed is the same as in a Brinkmann fluid (black solid 
curve) for solvent-mediated network-sheet BC when the net- 
work does not deform (small rjuj/fi). Blue and magenta curves 
are for direct network-sheet BC, H = 0, and A = 0; red and 
green curves are for solvent-mediated network-sheet BC and 
A = 0. Red (long dashes): rju/^ — 0.01. Green (medium 
dashes): rju/fi — 10. Blue (short dashes): rjuj/^i = 0.01. Ma- 
genta (dash-dots): 'qu:/^ — 10. 



modulus is « 4 Pa and viscosity is « 20 Pa s [27]. Esti- 
mating r w Jy/C^; where ^ is the mesh size, and using a 
wavelength of sperm flagellar beating (w 10 /xm) to deter- 
mine /c, we estimate that near ovulation rjuj/ ^ ss 10 and 
T/{r]k^) sa 0.03, while away from ovulation rjcu/ij, w 100 
and T /{rfk'^) w 10. Although for sperm the swimming 
speed enhancement seen in the stiff network/low frequency 
regime is unlikely to be realized in cervical mucus, for 
other swimmers in gels it may still be an important effect. 

As we were completing this letter, we learned of the 
work of H. Wada [15], who has independently studied a 
closely related problem. 
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